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ABSTRACT 

We present an analysis of the optical spectra of a volume-limited sample of 375 
radio galaxies at redshift 0.4 < z < 0.7 from the 2dF-SDSS Luminous Red Galaxy and 
QSO (2SLAQ) redshift survey. We investigate the evolution of the stellar populations 
and emission-line properties of these galaxies. By constructing composite spectra and 
comparing with a matched sample of radio-quiet sources from the same survey, we 
also investigate the effect on the galaxy of the presence of an active nucleus. 

The composite spectra, binned by redshift and radio luminosity, all require two 
components to describe them, which we interpret as an old and a younger population. 
We found no evolution with redshift of the age of the younger population in radio 
galaxies, nor were they different from the radio-quiet comparison sample. Similarly, 
there is no correlation with radio power, with the exception that the most powerful 
radio sources (Pi. 4 > 10 26 W/Hz) have younger stars and stronger emission lines than 
the less powerful sources. This suggests that we have located the threshold in radio 
power where strong emission lines "switch on", at radio powers of around 10 26 W/Hz. 
Except for the very powerful radio galaxies, the presence of a currently-active ra- 
dio AGN does not appear to be correlated with any change in the observed stellar 
population of a luminous red galaxy at z ~ 0.5. 
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1 INTRODUCTION 

It is now accepted that most, if not all, massive galaxies 
host super-massive black holes at their centre. The observed 
tight correlation between the mass of the black hole and the 
global properties of t he host galaxy l|Magorrian et al1ll998l ; 
iGebhardt et~ai]|2000l ) suggests that the growth of the black 
holes at their centres is intimately related to the formation 
and assembly of these galaxies. However, at any one time 
only ~5% of galaxies are seen as powerful radio sources or 
show other evidence of an active galactic nucleus (AGN). 
The question arises: are the host galaxies of such AGN dif- 
ferent from the hosts of inactive sources? 

It has been known since the early 1960s that the 
hosts for powerful radio galaxies are massive ellipticals 
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(e.g. iMatthews et all 1 1964 ). The trigger for the transition 
from quiescence to an active state has been suggested to 
be gravitational interaction between galaxies. Signatures of 
tidal interactions are common amongst radio galaxies (e.g. 
iHeckman et al.lll986l h which presumably result in a large 
increase in the amount of material feeding the black hole, 
thus triggering the radio emission. Some individual radio 
galaxies also show s ignificant association with young stel- 
lar populations (e.g. lAretxaga et alj 200 ll: Tadhunter et al.1 
120021 : 1 Wills et al.l [20021 : 1 Johnston et al.ll2005l h this suggests 
that the merger process which triggers the radio emission 
can also lead to star formation in the host galaxy. However, 
numerical studies of galaxy interactions and mergers show 
that star formation is not necessarily enhanced , due to the 
amount of gas swept out (di Ma tteo et al.| [2007). 

Several questions remain to be answered in this sce- 
nario: What is the relationship among mergers, the trigger- 
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ing of the active nucleus, and star formation: are all three 
components necessarily present? Do all early-type galaxies 
go through radio galaxy phases as they assemble? 

In order to investigate these questions, we need to study 
a large and uniform sample of radio galaxies. Samples of 
this kind h ave recently been stud ied in the local univers e 
(z < 0.3) bv lBest et all (|2005h and lMauch fc Sadled i|2007f) . 
but until now there has been no comparable sample available 
at hi gher redsh ift. 

iBest et all (|2005f ) investigated the properties of a sam- 
ple of 2215 nearby radio galaxies (most with 1.4 GHz radio 
luminosity below 10 25 W/Hz) from the Sloan Digital Sky 
Survey. They showed that the fraction of galaxies hosting 
a radio-loud AGN was a strong function of stellar mass, 
and that there was no correlation between radio luminosity 
and opt ical emission-line luminosity for the galaxies in their 
sample. IBest et all l|2005h concluded that optical AGN and 
low-luminosity radio-loud AGN are independent phenom- 
ena which are triggere d by different physical mechanisms. 
iMauch fc Sadlerl i|2007h studied a sample of 2661 radio-loud 
AGN selected from the 6dF Gal axy Survey. They confirmed 
the findings of IBest et ail (|2005l ) that radio-loud AGN pref- 
erentially inhabit the brightest and most massive host galax- 
ies, and showed that the fraction of all galaxies which host a 
radio-loud AGN scales with the infrared K-band luminosity 
as . 

We have constructed a volume-limited sample of radio 
sources by combining optic al data from the 2SLAQ LRG 
survey (| Cannon et a l. 2006) with the r adio data from the 
VLA FIRST survey (|Becker et al.lll995l ). The 2SLAQ sam- 
ple is focused on luminous red galaxies, thus pre-selecting 
for the large elliptical galaxies m ost likely to host powerful 
radio sources. ISadler et al.l (|2007t ) used this sample to study 
the evolution of the radio sources; in the current work, we 
investigate the properties of the host galaxies. 

The signal-to-noise ratio of individual 2SLAQ spectra 
is adequate for redshift determinations, but is generally too 
low for accurate measurement of weak spectral features. We 
overcome this limitation by combining spectra to form com- 
posites. This approach allows us to investigate the variation 
of composites with both redshift and with radio power, and 
to address the key question: How does the stellar popula- 
tion of a galaxy correlate with the presence or absence of an 
active radio-loud AGN in its centre? 

Throughout this paper, we use Ho = 71 km s~ Mpc - , 
Q. m = 0.27 and Q. A = 0.73. 



2 DATA REDUCTION 



2.1 Source selection 



The source selection is discussed in detail in ISadler et al.l 
l|2007i . Paper I); we present a brief summary here. 

The 2SLAQ Luminous Red Galaxy survey 
( Oannon et all 120061 ) consists of 14,978 spectra taken with 
the Two-degree Field instrument (2dF) on the 3.9m Anglo- 
Australian Telescope , of sources selected to have ugriz 
|Fukugita et al.l|l996l ) photometry a nd non-stellar sha pe in 
the Sloan Digital Sky Survey (SDSS: lYork et al.ll2000h . The 
catalogue of sources with spectr a was then cross-matched 
with the VLA FIRST catalogue l|Becker. White fc Helfandl 



Il995l ) and the NVSS (|Condon et al.lll998h . These two sur- 
veys have complementary properties, with NVSS sampling 
the total flux density of extended sources, while FIRST 
has higher spatial resolution. Galaxies within 30 arcsec of 
a FIRST/NVSS source were flagged as possible matches 
(giving 2871 potential candidates). Optical images of the 
galaxies, p rimarily from Sloan D igital Sky Survey Data 
Release 3 (jAbazaiian et all [2005), with a small number 
taken from the SuperCOSMOS Sky Surveys l|Hamblv et al.l 
2001), were compared visually with radio images from the 
FIRST survey. This process yielded 391 optical galaxies 
with 2SLAQ spectra which are true associations with radio 
sources. During this visual inspection, we also flagged and 
accounted for sources which are double and/or complex in 
structure. 

Of these 391 radio sources, 14 had redshift quality val- 
ues Q < 3, indicating a less-reliable redshift determination. 
These objects were excluded from further analysis; the re- 
mainder of the objects, with Q > 3, have redshift reliability 
> 95% (see lCannon et alfeOOrl . §5.4). The total radio detec- 
tion rate was 2.7%. There were two sources in this list whose 
spectrum showed a measurable redshift but was significantly 
contaminated by a foreground sta r (the sources indic ated as 
"M star plus galaxy" in Table 3 of lSadler et"al]|2007l ). While 
these were usable for constructing a radio luminosity func- 
tion, they would complicate our spectral modelling, so were 
removed from the sample for this paper, leaving 375 radio 
sources with optical spectra. 

The total NVSS/FIRST radio flux of each object, in- 
cluding the flux from all separate components identified by 
eye, was converted to total radio power assuming Ho = 
71 kms _1 Mpc _1 and Qm = 0.27, using a fe-correction of 
the form (1 + z)~^ 1+ a ' with a = —0.7, w hich is the median 
radio spectral index l|Mauch et al.|[2003l) . 

The final source list of 375 objects has redshifts ranging 
from z = 0.31 to 0.76, and radio power between Pi.4GHz = 
1Q 23.44 Wj/Rz and 1Q 27.02 W / Hz The redshi ft distribution 

of radio sources is shown in Fig. [T] The median redshift is 
z = 0.54, and the median radio power is logP = 24.60. 

2.2 Creation of a matched sample 

To investigate the effect which the presence of an active 
nucleus has on its host galaxy, we must create a sample 
of non-radio galaxies, matched as closely as possible to the 
radio-loud sample. Since the probability that an early-type 
galaxy is a radio source is a strong function of op tical lumi- 
nosity (|Auriemma et alH "l977: Sa dler et al.ll 1989] ), the radio 
galaxies in our sample are significantly more luminous than 
the 2SLAQ LRG sample as a whole, which makes careful 
matching imperative. 

For each radio source, we selected a matching source 
from the remainder of the 2SLAQ catalogue. The sensitivity 
of the FIRST ra dio catalogue (with a typical flux density 
limit of 0.75 mjy; iBecker et aTlll995l ) means that the upper 
limit to the radio power of these "radio-quiet" sources range 
between 10 23 and 10 24 W/Hz at 1.4 GHz. For each radio- 
loud object, we found the object best matched in redshift 
and r-magnitude by minimising the quantity 



d = (zrl 



zhq) 2 + (r RL 



r RQ ) 2 /20 



(1) 



There are several possible ways of finding the "nearest" 
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radio-quiet galaxy; this expression was chosen empirically to 
weight the magnitude and redshift equally, so the resulting 
sample has roughly the same distribution in redshift and 
apparent magnitude. Since radio-loud sources make up only 
2.7% of the LRG sample, this approach produces a sample 
which is well-matched to the radio-loud sources in redshift 
and apparent r magnitude. The maximum discrepancy be- 
tween a radio source and its "matched" galaxy was 0.18 mag 
in apparent magnitude and 0.0067 in z, with 90% of sources 
having Ar <, 0.03 mag and Az <, 0.001. 



2.3 Spectral reduction 

We began with the wavelength-calibrated spectra produced 
by the 2SLAQ project. These had a typical wavelength cov- 
erage of 5000-7250A, which means that, at the redshift 
range of our targets, the region containing Call H and K 
and the 4000A break is always covered. Although our spec- 
tra go well below 4000 A rest wavelength, our observational 
data are in the spectral region where instrumental and at- 
mospheric effects on continuum slope are not severe. For 
redshifts beyond z ~ 0.45, H/3 and [Om] 5007 are shifted 
out of the 2dF wavelength range, so the most prominent 
emission line is [O n] 3727. 

To combine the spectra, several steps of processing were 
applied to each spectrum. The regions around the strongest 
night sky features were removed and replaced by an interpo- 
lated value. We are hampered by the fact that the 2SLAQ 
spectra are not flux-calibrated, which means that the instru- 
ment response can introduce changes in the slope of the con- 
tinuum. We attempted to correct for this, at least in a mean 
sense, by dividin g by an average re s ponse curve. We use the 
curve derived by iRoseboom et al.l (|2007T ) , which compared 
the 2SLAQ spectra of 160 galaxies which were also observed 
by the SDSS LRG survey (Eisenstein et al. 2001, 2003). 
We divided each spectrum by this response curve: since we 
are producing composite spectra this mean response curve 
should suffice, even though it does not take into account 
night-to-night variations. We then shifted each spectrum to 
the rest frame using the known redshift, and rebinned all 
the spectra onto the same scale, between 3280-4870A, with 
a dispersion of 1.4A/pixel 



2.4 Creation of composites 

The distribution of sources in (z,logP) is shown in Fig. [T] 
Although the redshifts range from z — 0.31 to 0.76, and the 
radio power from logP = 23.44 to 27.02, the bulk of sources 
cluster near the middle of the range in each parameter. We 
thus had two choices: we can create composites with equal 
numbers of sources, which produces composites with equal 
signal-to-noise, or we can create composites within equal 
intervals in redshift or radio power, which maximises the 
baseline over which we can investigate changes, but at the 
expense of lowering the signal-to-noise ratio of the outlying 
composites. 

We thus created four separate sets of composite spectra 
(Table [T] shown graphically by the dashed lines in Fig. [TJ : 

(A) sort spectra by redshift and create 5 composites with 
equal numbers of sources: each composite is composed of 75 
sources, and the redshift limits are shown in Table [1] 
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Figure 1. Distribution of sources in redshift and radio power at 
1.4 GHz (measured in W/Hz). The separate histograms of redshift 
and radio power are shown above and to the right, respectively. 
The binning for the four separate sets of composites are indicated 
by the lines. Set A, binning by redshift with 75 sources in each 
bin, is shown by the vertical blue dashed lines; set B, binning by 
redshift with equal spacing, is shown by the vertical red dotted 
lines. Set C, binning by radio power with 75 sources in each bin, 
is shown by the horizontal blue dashed lines; set D, with equally 
spaced bins in radio power, is shown by the horizontal red dotted 
lines. 



(B) sort by redshift and create 4 composites with equal 
spacing in redshift: this gave us four groups, with z < 0.45 
(28 sources), 0.45 < z < 0.55 (196 sources), 0.55 < z < 0.65 
(138 sources), and z > 0.65 (28 sources); 

(C) sort by radio power and create 5 composites with 
equal numbers of sources: each composite is composed of 75 
sources, and the redshift limits are shown in Table [T] 

(D) sort by radio power and create 4 composites with 
equal spacing in radio power: this gives us four bins, with 
logP < 24 (27 sources), 24 < logP < 25 (269 sources), 
25 < logP < 26 (83 sources), and log P > 26 (11 sources). 

There are several different approaches which may be 
taken in combining low signal-to-noise ratio spectra to form 
composites: sp ectra may be weighted by the variance in e ach 
spectrum (e.g. lFrancis et aljl9 91; Eisenstein et al. 20 Q3|), or 
weigh ted by luminosity ( Baldrv et al.ll2002l : iDressler et al.l 
2004); or the composite can be constructed from the me- 
dian of individual spectra at each wavelength l|Croom et al.l 
2002). Since the 2SLAQ galaxies are all intrinsically lumi- 
nous red galaxies, we have taken the simplest approach and 
created composites by averaging the spectra in each group. 
The spectra were all taken with very similar wavelength cov- 
erage but have different redshifts, so in each rest wavelength 
bin, different numbers of spectra contribute to the sum, par- 
ticularly for the highest and lowest redshifts. We corrected 
for this by counting the number of spectra which contribute 
at each wavelength and dividing the composite by this num- 
ber. 

The composite spectra were then normalised by the flux 
at 4050A. Since the resolution of the spectral models we 
chose is 10A (see § I3.2[l , we binned each composite into 20A 
bins, and then removed bins which included the bright emis- 
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Table 1. Parameters of the composites and results of the fits. Four sets of composites A-D were created (see £j I2,4l for details). Columns 
f-4 show the name of the composite set, the range in redshift or radio power, and the number of spectra which went into each composite. 
The next three columns show the results of the spectra modelling (§ 13. 2^ : the age of the young population, the fraction of light contributed 
by the young population at 4050A, and the reduced xl value for the best fit (for v = 31 degrees of freedom). Corresponding composites 
of a matched sample of radio-quiet sources (§ 12.21 1 were created, and the results of fitting to these are shown in columns 8-10. 

Radio-loud sample Radio-quiet sample 

No. of Age Age 



Sample 


# 


Range 


galaxies 


(Myr) 


Fraction 


xl 


(Myr) 


Fraction 


xl 




1 


2 = 0.31-0.48 


75 


690 


0.40 


3.55 


710 


0.54 


4.42 


(A) Redshift: 


2 


2 = 0.48-0.52 


75 


1000 


0.47 


3.61 


1130 


0.56 


2.48 


equal numbers 


3 


2 = 0.52-0.56 


75 


700 


0.45 


3.29 


690 


0.46 


1.84 




4 


2 = 0.56-0.59 


75 


1120 


0.55 


1.32 


670 


0.50 


0.87 




5 


2 = 0.59-0.76 


75 


1130 


0.60 


1.32 


1120 


0.70 


1.55 




1 


2 = 0.31-0.45 


26 


690 


0.35 


6.48 


1110 


0.57 


2.00 


(B) Redshift: 


2 


2 = 0.45-0.55 


192 


700 


0.40 


3.26 


1000 


0.54 


3.19 


equal spacing 


3 


2 = 0.55-0.65 


133 


1120 


0.66 


1.26 


1120 


0.66 


1.52 




4 


2 = 0.65-0.76 


24 


1130 


0.64 


0.94 


1130 


0.81 


0.55 




1 


logP = 23.44-24.29 


75 


720 


0.50 


3.61 


710 


0.53 


2.68 


(C) Radio power: 


2 


logP = 24.29-24.51 


75 


1120 


0.52 


1.90 


1010 


0.59 


1.42 


equal numbers 


3 


logP = 24.51-24.75 


75 


710 


0.45 


1.94 


710 


0.64 


6.45 




4 


logP = 24.75-25.09 


75 


700 


0.47 


2.06 


1130 


0.58 


0.90 




5 


logP = 25.09-27.02 


75 


1140 


0.53 


1.06 


700 


0.56 


1.90 




1 


logP = 23.44-24 


26 


710 


0.55 


3.81 


710 


0.45 


3.00 


(D) Radio power: 


2 


log P = 24-25 


259 


700 


0.46 


2.84 


710 


0.55 


2.84 


equal spacing 


3 


log P = 25-26 


79 


1130 


0.51 


0.81 


990 


0.64 


1.77 




4 


logP = 26-27.02 


11 


110 


0.46 


1.61 


680 


0.59 


0.68 



sion lines, principally [On] and [Om]. We estimated the er- 
rors on the flux in each wavelength bin by measuring the 
RMS of the 14 original pixels which were included in each 
20 A bin. We chose to perform the fits to the models over 
the wavelength range 3500-4400A, based on our desire to 
fit the same wavelength region in all spectra; a total of 33 
independent data points in each spectrum were used in the 
fits. 

Composites for the matched sample of radio-quiet 
sources (§ I2.2p were created in exactly the same manner; 
this data set serves as a control sample to see how the stellar 
population varies in the absence of a powerful radio source. 



3 DATA ANALYSIS 

3.1 Differences between composites 

To determine if the presence of an active nucleus affects the 
galaxy, we compared the composites formed from the radio- 
loud sample and the matched radio-quiet sources (§ 12.21) . 
The most basic test is to use a simple \ test to find the 
probability that a pair of composites come from the same un- 
derlying distribution. Each of the eighteen radio- loud com- 
posite spectra (Table [IJ was compared with the correspond- 
ing radio-quiet counterpart. All composites were consistent 
with being from the same distribution, except the D4 com- 
posite formed from the highest-power radio sources, where 
the xl — 3.5. This initial evidence suggests that, with the 
exception of the very highest-power radio sources, there is no 
evidence of difference between galaxies containing an active 
AGN and those that do not. 



3.2 Modelling the stellar continuum 

We modelled the observed continuum as an old population 
plus a blue component due to a second population of stars 
of a younger age, using mo dels from the GISSEL96 library 
|Bruzual A. fc Charlotlll993r ). Of the models in their library, 
we used the ones with a Salpeter initial mass function with 
ma ss limits of 0.1 and 100 Mq and solar metallicity, using 
the lGunn fc Strvkerl l|l983t ) stellar spectral atlas. 

We used two single-age stellar populations, representing 
instantaneous bursts of star formation. The old population 
was represented by a single-age population of age 7000 Myr, 
the younger stars by a population with a range of ages (10, 
20, 50, 70, 100, 200, 500, 700, 1000, 2000 or 5000 Myr). The 
exact age of the "old" population is not critical, since the 
spectrum changes little at these ages; we chose 7 Gyr since 
at the highest redshift in our sample [z = 0.76) the age of 
the universe was 7 Gyr. We assembled the model spectra 
from the GISSEL96 library, and binned them to the same 
resolution as the observed spectra. 

Since the observed continuum is modelled as the sum of 
two populations, we assumed each composite is represented 
as a 7000 Myr population plus some fraction / of a younger 
population with age r; our task is to determine / and r 
for each composite. This parametrisation is unlikely to be 
a realistic representation of the stellar populations of these 
galaxies. However, the exact form of the models is not im- 
portant, because our study is a differential test, comparing 
the properties of the radio-loud galaxies with their radio- 
quiet counterparts. 

We created a 11 x 11 grid of model spectra. Along one 
axis, the age of the young population r is varied between 
10 Myr and 5 Gyr; along the other axis the fraction of light 
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A composites B composites 




0.4 0.5 0.6 0.7 0.4 0.5 0.6 0.7 

z z 

Figure 2. Results of fits to the A and B composites, showing 
the age and fraction of light contributed by young stars as a 
function of rcdshift. The red filled symbols (solid lines) represent 
the radio-loud sources, while the blue open symbols (dashed lines) 
represent the matched sample of radio-quiet sources; the radio- 
quiet sources have been offset slightly in redshift for clarity. The 
age of the young population does not change with redshift; the 
fraction of light contributed rises slightly. 



contributed by the young population varies linearly between 
and 100%. Thus the spectrum at grid position (5,3) con- 
sists of an old (7 Gyr) population plus a 100 Myr population, 
where the young population is contributing 20% of the light 
at 4050 A. 

We calculated the x 2 difference between the model spec- 
trum and each of the observed composite spectra for all 121 
model spectra. We allowed the normalisation of the model 
spectrum to vary between 0.75 and 1.25 of the flux of the ob- 
served spectrum at 4050A, to allow for possible differences 
in slope between the model and the observed spectrum. This 
yields a \ 2 map for each observed composite; the best fit is 
chosen to be that where the x 2 value is minimised. 

The results of the fitting are shown in Table Q] and 
shown graphically in Figures [2] and [3] The best-fitting age 
and fraction were derived from a parabolic fit to the one- 
dimensional x 2 curves through the minimum. The error bars 
represent the confidence interval given by the ellipse repre- 
senting Ax 2 = 3 (containing ~ 90% of normally distributed 
data). In all cases, the x 2 OI the two-component model fit 
was significantly (~ factor 2) better than the x 2 °f the fit 
using a single-aged stellar population. 

The results of our modelling can be summarised as fol- 
lows. There is no evolution with redshift of the stellar popu- 
lation of the radio galaxies, nor are they in any way different 
from radio-quiet galaxies matched in redshift and r luminos- 
ity. The age of the "young" population is unchanged with 
redshift, with our spectra requiring a second population of 
stars of age ~ 820 Myr over the whole redshift range of our 
sample. The fraction of light contributed by young stars is 
about 50%, and while it rises slightly over the redshift range 
of the sample, the change is not significant (1.4cr). 

The B-composites and their best-fitting models are 
shown in Fig. [4] 

We can also examine how the spectrum changes as a 
function of radio power, using the C and D composites de- 
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Figure 3. Same as Fig. [2] but for the C and D composites, show- 
ing the age and fraction of light contributed by young stars as a 
function of radio power. Neither the age of the young population 
nor the fraction of light contributed changes significantly as a 
function of radio power, except at the highest radio power (com- 
posite D4), where the second population is significantly younger. 
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Rest wavelength (A) 

Figure 4. Model fits to the four B composite spectra, with equal 
spacing in redshift. The composite spectrum is shown in red, with 
the model fit, consisting of an old population plus a young pop- 
ulation of age t contributing a fraction / of the light at 4050A, 
shown as the thick black line (Table [TJ. The dotted lines indicate 
the wavelength range over which the fitting was performed. 



fined earlier. Again, our modelling shows no change in the 
age of the population or the fraction of young stars as the 
radio power increases (Fig. [3}. The only exception is that 
the very highest power sources (the D4 composite, consist- 
ing of the 11 sources with logP > 26; Table [l} appear to 
have a substantially younger population (100 Myr), as well 
as strong [O n] 3727 emission (which was not included in 
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Figure 5. Model fits to the four D composite spectra, with equal 
spacing in radio power. The composite spectrum is shown in red, 
with the model fit shown as the thick black line (Table [TJ- Note 
the substantially stronger [O n] emission line in the highest radio 
power spectrum (composite D4), as well as [Nem] emission. 



the fit). This difference is not seen in the matching sam- 
ple of radio-quiet sources, suggesting that the young stars 
are in some way associated with the active nucleus. Recall 
that this was the only composite showing a simple \ 2 dif- 
ference with the other composites (§ I3.1|l . We checked the 
eleven individual spectra which went into the D4 composite, 
in case the composite was being skewed by a single peculiar 
object. No individual spectrum appeared to be peculiar; 6 
of the 11 showed strong emission lines. The D-composites 
and their best-fitting models are shown in Fig. [S] We tried 
splitting the second-highest power bin (D3), to see if there 
was a continuous trend; both sub-bins had similar ages, of 
~ 1000 and ~ 700 Myr, quite different from the 100 Myr 
population found in the D4 bin. Thus the difference appears 
to be confined to the very highest radio-power sources. 

In Figure [6] we plot the ratio of these models. The fig- 
ure shows there are small changes to spectral lines, but the 
principal difference is a change in slope in the model fit to 
the D4 composite. 

We investigated whether the young population in the 
highest-power radio sources could in fact be due to emission 
from the central AGN, which would also produce excess blue 
light. We tried replacing the second, young population in 
our fits with a power-law continuum of the form f\ oc X a , 
where a was allowed to vary between —5 and 0. The fit was 
significantly poorer than the best-fitting model with a young 
stellar population: the 4000A break was not well reproduced. 
Thus we conclude that the highest power radio sources do 
show evidence of association with a population of young 
stars. 
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Figure 6. Comparison between the model fits to the four D com- 
posite spectra (Fig. [B). The models are shown divided by the 
best-fit model to the Dl composite (r = 500 Myr, / = 0.5), and 
offset in the vertical direction for clarity. 

3.3 Emission lines 

Emission lines were present in many of the spectra, mostly 
[O ill, and occasio nally [Nem] and [Nev] (see Table 3 of 
ISadler et all \200H) . The [O ill] 4958,5007 pair was seen 
in only a single galaxy, J100322.41-000137.8, but was 
redshifte d beyond th e wave length range of most of our 
sources. ISadler et al l (|2007h found that the fraction of 
2SLAQ radio galaxies which show [O n] emission is higher 
(27%) than the overal l 2SLAQ spectroscopic sample (17.7%; 
iRoseboom et al . 2006). However, since more luminous galax- 
ies are more likely to show emission, this result needs to be 
approached with caution. In any case, most of the radio 
galaxies in 2SLAQ would not have been recognised as AGN 
on the basis of their optical spectra alone. This agrees well 
with several other studies, which found little evidence of cor- 
relation between radi o power and emiss i on line strength i n 
nearby radio galaxies l|Rixon et al.lll99l1 ; IOwen et alJll995h . 

By using a comparison sample of radio-quiet galaxies, 
matched in redshift and luminosity, we are able to investi- 
gate this effect. Comparison of the emission- line properties 
of our radio-loud sources with the matched sample of radio- 
quiet galaxies reveals differences between the two groups. 
We can examine the emission lines in two ways: finding the 
number of individual galaxies which show emission lines as a 
function of various parameters, and investigating the prop- 
erties of the emission lines in the composite spectra. 

3.3.1 Emission lines in individual sources 

The equivalent widths of the [O n] line in the individual spec- 
tra were measured using the SPECFIT package implemented 
in I RAF (Kriss 1994). We normalised the spectra by a spline 
fit to the continuum, then measured the equivalent width of 
the emission line by fitting a Gaussian profile to the line. We 
must be careful to account for sources where the emission 
line fell on a sky line, or was redshifted out of the spectrum, 
and so would not have been detected if it were present. There 
were 333 sources in which [O n] was potentially detectable, 
with sources at redshifts ~ 0.48 and ~ 0.58 most strongly 
affected when the [O h] line fell on a bright sky line. 

Following ISadler et all l|2007h . we used a rest equiva- 
lent width cutoff of 7A to define sources with [O n] emis- 



© 0000 RAS, MNRAS 000, 000-000 



Radio galaxies in 2SLAQ 7 



sion. Using this criterion, 50 out of the 333 (15%) radio-loud 
sources in which [O n] emission was potentially detectable 
show such emission, compared to 25 out of the matched 
radio-quiet sample (7.5%); so the radio-loud sources are sig- 
nificantly more likely to show [O n] emission than galaxies 
matched in optical luminosity and redshift. If we assume the 
fraction of radio-quiet galaxies showing emission represents 
the null-hypothesis (25/333), then observing 50/333 radio- 
loud sources with emission is highly significant (p <C 0.001). 
For the [Nem] 3869 and [Nev] 3426 lines, we used a cut- 
off of W\ > 4 A to define the emission-line sources. There 
were 25 sources with [Nem] emission and 15 with [Nev] 
emission in the radio-loud sample, compared to 11 and 8 
respectively in the radio-quiet sample. Again, the radio- 
loud sources are roughly twice as likely to show emission 
than their radio-quiet counterparts (probability p < 0.001 
and 0.05 respectively). Twenty-one radio-loud sources have 
IF\([0 m]) > 15A, compared with just five sources among 
the radio-quiet galaxies. 

There is a slight trend for the fraction of sources show- 
ing emission to increase with redshift. In Fig. [7] we plot 
the fraction of sources which show [0 11] emission (out of 
those in which it was potentially detectable) as a function of 
redshift. A similar trend is seen in the radio-quiet sources. 
This could be an absolute magnitude effect, as we are se- 
lecting slightly mo re luminous galaxies at higher redshift 
(|Sadler et al.ll2007l ). 

[0 11] emission is seen in radio galaxies at all radio pow- 
ers (though it is somewhat more common in low- (logP < 
24.5) and high-power (log P > 25) sources, as shown in 
the right-hand panel of Fig. [7}. "Strong" [On] emission, 
however, with W\ > 15A, is significantly more common 
in the highest-power radio galaxies: 28% of sources with 
log P > 25.5 (9 out of 35) have W\ > 15 k compared to 
~ 3% of sources at lower radio power. [Nem] and [Nev] 
emission is also more common at higher radio power. Nei- 
ther of these trends can be seen in the sample of radio-quiet 
sources which match the distribution of the radio galaxies 
in z and r. 

3.3.2 Emission lines in the composite spectra 

By combining many spectra, we can measure changes in 
line ratios much more accurately than we can in individual 
sources. 

We measured the line properties of the emission lines 
in the composite spectra by fitting Gaussian profiles using 
SPECFIT again, using the same technique as in t|3.3.1l The 
spectra were normalised by a spline fit to the continuum, and 
then Gaussians were fit, constraining the wavelengths to be 
within ±5A of the known wavelengths of the lines. These 
fits are shown in Figure [8] the line properties are listed in 
Tabled 

There is little difference in the emission lines of the com- 
posite spectra until the very highest radio-power is reached. 
This composite shows much stronger [O n] emission, as well 
as strong [Nev] and [Nem] emission lines, lines which are 
characteristic of quasars. The radio-quiet galaxies which are 
matched in luminosity and redshift show no increase in 
[0 11] emission, and no trace of neon lines. This is not just 
due to the contribution from a single source: three of the 
nine sources with highest radio power have VF3426 > 4A (in 
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Figure 7. Fraction of sources with [On] and [Nev] emission, 
as a function of redshift and radio power; log P is defined for 
the radio-quiet sample by the matching to the radio-loud sources 
(ii 12. 21 . The red filled symbols represent the radio- loud sources 
(RL), while the blue open symbols represent the matched sample 
of radio-quiet sources (RQ). The squares and solid lines show the 
fraction of sources with EW([Oll] 3727) > 7A, the triangles and 
dashed lines show the fraction of sources with EW([Ne v] 3426) > 
4A. The radio-loud sources show a higher fraction of sources with 
emission than the radio-quiet sources, as well as a tendency for 
this fraction to increase with redshift. The likelihood of showing 
emission increases with radio-power (right-hand panel), especially 
for the [Ne v] emission. 



Table 2. Equivalent widths of the fits to the emission lines in 
the D composites (Fig. [SJ The equivalent widths of the three 
emission lines are shown, for the radio-loud and matched radio- 
quiet composites. The last two lines show the fits to the Neon 
composite and the [O II] composite ( §3.3.21 . 





[On] 3727 


[Nev] 3426 


[Nem] 3869 


Composite 


(A) 


(A) 


(A) 


Radio loud Dl 


4.2 ±0.7 


< 1.0 


0.9±0.7 


Radio loud D2 


3.7 ±0.4 


1.0±0.4 


1.0 ±0.4 


Radio loud D3 


3.3 ±0.4 


< 0.5 


1.4 ±0.4 


Radio loud D4 


14.1 ±0.6 


6.8 ±0.5 


5.2±0.4 


Radio quiet Dl 


1.5 ±0.4 


< 0.5 


< 0.6 


Radio quiet D2 


2.1±0.3 


< 0.4 


1.2 ±0.4 


Radio quiet D3 


2.2 ±0.3 


< 0.3 


1.5 ±0.4 


Radio quiet D4 


1.8±0.5 


1.9 ±0.9 


1.0 ±0.9 


Neon composite 


17.2 ±0.7 


7.9 ±0.6 


8.1 ±0.6 


[O II] composite 


17.7 ±0.9 


< 0.7 


1.8 ±0.6 



two sources the line fell on a sky line). This result implies 
that the increased emission seen in the most powerful radio 
sources is not just a function of galaxy size, but is tied to 

the activity of the cent ral black hole. 

lOwen et all (119951 ) and lBest et~ai1 i|2005h . found no cor- 
relation between radio power and emission line properties; 
their samples only reached radio powers P ~ 10 25 W/Hz. 
The emission-line luminosity of very power ful radio sources 
is known to correlate with radio power (e.g. lMcCarthvlll993l , 
Fig. 2). There must be a threshold radio power at which 
strong emission lines "switch on" at a level above the usual 
incidence in quiescent galaxies; at the highest radio lumi- 
nositi es, virtually every rad io galaxy has strong forbidden 
lines (|Hine fc Longairlll979h . Our results suggest we have 
found this threshold, at radio powers of around 10 26 W/Hz. 
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01: log P = 23.63-24 



D2: log P = 24-25 



Aa 



D3: log P = 25-26 



D4: log P = 26-27.02 
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Figure 8. Gaussian fits to the emission lines in the D compos- 
ites, with equal spacing in radio power. The normalised composite 
spectra are shown in red, and the Gaussian fits to the emission 
lines are shown as the thick black line. [O o] emission is seen at 
all radio powers, though it is significantly stronger in the highest 
power sources; these are the only ones to show significant [Nov] 
and [Ne Hi] emission. The fourth line is a weaker [Ne v] line at 
3346A. The blue end of the lowest-power composite is noisy, be- 
cause the lowest-power sources are only found at nearby redshifts, 
so there is not much signal at short wavelengths in this composite. 



As noted earlier (§ 13.3.10 . there were 21 sources which 
showed strong [O n] emission, with W\ > 15A. There were 
also 32 sources which showed [Nem] and/or [Nev] emis- 
sion. There was not total overlap between these: there 
were ten sources which showed strong [O n] emission but 
no [Ne m] /[Nev] emission. To test whether this difference 
was real, we constructed two composite spectra: (a) the 32 
sources which showed [Nev] and/or [Nem] emission (the 
"Neon composite"), and (b) the sources that showed strong 
[On] emission but no [Nem] or [Nev] emission (the "[On] 
composite"). We could only use five of the ten sources with- 
out [Ne ill] /[Nev] emission in constructing the [On] com- 
posite, because the other five had at least one of the three 
neon lines falling on a sky line. The composites are shown 
in Figure [10] 

Clearly these two composites have significantly different 
line ratios: the non-detection of [Ne v] in some sources is not 
just a signal-to-noise ratio issue. The Neon composite has 
line ratios [Nev] /[On] = 0.46 ± 0.03 and [Nem]/[0 n] = 
0.47 ± 0.03, while the [On] composite has [Nev] /[On] < 
0.04 and [Nem]/[0 n] = 0.10 ± 0.03. This suggests different 
ionisation mechanisms are operating in these two groups. For 
comparison, the composite quasar spectrum published by 
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Figure 9. Same as Fig. [8] but for the matched radio-quiet com- 
posites. Clearly the emission line properties are very different to 
the radio-loud sources, with no detectable [Nev] or [Nem] emis- 
sion. 
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Figure 10. Gaussian fits to the emission lines in composites of 
sources which show significant [Nev] and [Nem] emission, and 
sources which show strong [O n] emission but no [Ne v] or [Ne III] 
emission. Only sources where all emission lines were clear of the 
night sky lines were included. 



IVanden Berk et all (|200lh has ratios [Nev]/ [On] = 0.78 ± 
0.02 and [Nein]/[0 n] = 0.88 ± 0.02. 

We investigated whether the stellar continuum in these 
emission line composites were at all different; nothing defini- 
tive could be found. The [O n] composite required a higher 
proportion of young stars to model the continuum than the 
Neon composite (80% compared to 40%), but this was only 



© 0000 RAS, MNRAS 000, 000-000 



Radio galaxies in 2SLAQ 9 



marginally significant; there were only five spectra in the 
[O n] composite. 



4 CONCLUSIONS 

We have investigated the optical properties of a volume- 
limited sample of radio galaxies in the 2SLAQ Luminous 
Red Galaxy survey. We created composite spectra in order 
to investigate how the stellar populations of these galaxies 
depend on redshift or radio power. By constructing a com- 
parison sample of galaxies matched in redshift and luminos- 
ity, we can eliminate the possibility that any effects we find 
are a result of the fact that radio-loud galaxies are on the 
average more luminous than radio-quiet galaxies. 

We modelled the underlying starlight using single-age 
stellar population models. None of the spectra were ade- 
quately described by a single model; they all demanded an 
additional stellar component. We achieved satisfactory fits 
assuming that each spectrum consisted of an "old" (7 Gyr) 
population plus a fraction of light contributed by a "young" 
population, whose age was allowed to vary between 10 Myr 
and 5 Gyr. Neither a power-law component nor a model with 
continuous star formation described the spectra as well as 
the two-stellar population model. 

The stellar populations of all the composites were al- 
most identical, independent of redshift, radio power, and 
whether or not the galaxy contains an active radio source. 
The fraction of light contributed by the young population 
(~ 40% at 4050A) increased marginally with redshift, but 
is not correlated with radio power. The age of the "young" 
component was consistent with an age of IOOO^jqq Myr 
for all composites. The sole exception was the composite 
made from the eleven most powerful radio sources, with 
log P > 26, which was fitted with a young population of age 
100 Myr. This was not seen in the corresponding composite 
of radio-quiet sources, which means it is unlikely to be due to 
a systematic effect, which would affect both radio-loud and 
radio-quiet galaxies equally. In addition. ISadler et alj (I2007T ) 
showed that the median redshift of 2SLAQ radio galaxies is 
independent of radio power, so there are unlikely to be sub- 
tle aperture effects affecting the sample. 

The two-population model is probably not a very real- 
istic description of the underlying stellar light. However, our 
study is a differential test, comparing the properties of the 
radio-loud galaxies with their radio-quiet counterparts. In a 
future paper, we will investigate other models to represent 
the stellar populations of these galaxies. 

The emission line properties of the population show 
more correlation with radio activity. Approximately 15% of 
the radio-loud sources show significant [O n] emission, com- 
pared to 7.5% of the radio-quiet sample; similarly, the radio- 
loud sources were roughly twice as likely to show [Nelll] 
and/or [Nev] emission than their radio-quiet counterparts 
(12% and 6%, respectively). Still, it is notable that most of 
the radio galaxies in the 2SLAQ would not have been recog- 
nised as AGN on the basis of their optical spectra alone. 
The different [Oll]/Ne line ratios seen in different objects 
suggest there is more than one emission mechanism operat- 
ing to ionise the gas: these are likely to include hot stars and 
post-AGB stars as well as ionisation by the central AGN. 

These observations support the model that all luminous 



red galaxies contain a central black hole, only some of which 
are active at any one time; the current activity of the central 
black hole, as measured by the radio emission, does not ap- 
pear to correlate with the stellar population in the galaxy. 
The sole exception appears to be that the most powerful ra- 
dio galaxies show more activity and younger stars than the 
rest of the population. Even so, with young (~ 100 Myr) 
stars contributing 30% of the light, the fraction of mass in 
young stars is only 0.5% of the mass of the galaxy. 

We need to be careful to distinguish between differences 
related to the activity of the central black hole and luminos- 
ity effects, since more massive galaxies are more likely to 
harbour AGN. We controlled for this bias by constructing 
a matched sample of radio-quiet galaxies. The comparison 
sample is matched in redshift and apparent magnitude, so 
the radio-quiet galaxies match the radio-loud sources closely 
in absolute luminosity. Since we saw differences in the most 
powerful radio galaxies, but not in their radio-quiet counter- 
parts, the difference most likely relates to the presence of an 
active black hole. In Paper I, the 2SLAQ radio galaxies were 
found to cluster m ore strongly than t he 2SLAQ LRG pop- 
ulation as a whole l|Sadler et al.ll2007l ). It was speculated in 
that paper that this might be a luminosity effect, since the 
2SLAQ radio galaxies are more luminous than the rest of 
the sample. However, recent work using samples matched in 
luminosity suggest this is not the explanation (Wake et al., 
in prep.). This suggests the radio galaxies really are prefer- 
entially located in denser en vironments , with more gas avail- 
able to feed the black hole. iBest et all l|2007h have recently 
suggested that a similar effect is seen in lower-redshift radio 
galaxies from SDSS. 

Our results support the findings of IBest et al l l|2005h . 
who found that the probability that a galaxy is radio-loud is 
independent of whether it is optically classified as an AGN, 
and that the host galaxies of radio-loud AGN match their 
radio-quiet counterparts. If optical emission lines trace cur- 
rent feeding of black holes, while radio emission points to 
mostly dormant black holes, accreting at low rates, then the 
current rate of growth of black holes in these galaxies is very 
low. 

One notable point is that the age of the "young" popu- 
lation remains unchanged at 700 Myr (Table [TJ despite the 
different ages of the galaxies: the highest redshift galaxies in 
our sample (z = 0.76) are 3.1 Gyr younger than the lowest 
redshift galaxies (z = 0.31). If the stellar populations in the 
lowest redshift galaxies were born 1 Gyr before z = 0.31, 
we should see the era of star birth, since it is covered in 
our sample (at a redshift z — 0.39). We see no evidence, 
however, for an epoch of star formation at that redshift. 

So where are the galaxies which are forming stars as 
we observe them? There are two possibilities. If the single- 
age stellar population is a bad assumption, we may have 
low-level ongoing star-formation instead of star formation 
in discrete bursts. The second possibility, however, is that 
the initial colour selection of the 2SLAQ sample is biased 
towards accepting galaxies with these ages; in other words, 
the counterparts of our z — 0.3 galaxies are indeed forming 
stars at z = 0.4, but they don't appear as red galaxies while 
forming stars, so are not included in the 2SLAQ sample. 

In fact, our modelling suggests that we may not notice 
these galaxies using the 2SLAQ gri colour selection criteria. 
A model galaxy with 5% of the mass in young stars formed in 
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a single burst (which contributes 40-50% of the light around 
1 Gyr after the burst) is only 1 mag bluer than the bulk of 
old galaxies. Furthermore, since most of the colour evolu- 
tion occurs very early on, in the first 50 Myr, the chances 
of finding such an object in our sample is not high. In or- 
der to locate these galaxies, we would need to assemble a 
large sample of luminous galaxies using a different colour 
selection. 



ACKNOWLEDGMENTS 

We would like to thank the referee, Jasper Wall, for his ex- 
tremely helpful comments and suggestions. Funding for the 
SDSS and SDSS-II has been provided by the Alfred P. Sloan 
Foundation, the Participating Institutions, the National Sci- 
ence Foundation, the U.S. Department of Energy, the Na- 
tional Aeronautics and Space Administration, the Japanese 
Monbukagakusho, the Max Planck Society, and the Higher 
Education Funding Council for England. The SDSS Web 
site is http : / /www . sdss . org/ . This work was partially sup- 
ported by National Science Foundation grant AST-0607634 
(N.P.R.). 



REFERENCES 

Abazajian K., et al. 2005, AJ, 129, 1755 

Aretxaga I., Terlevich E., Terlevich R. J., Cotter G., Diaz 

A. I., 2001, MNRAS, 325, 636 
Auriemma C, Perola G. C, Ekers R. D., Fanti R., Lari C, 

Jaffe W. J., Ulrich M. H., 1977, A&A, 57, 41 
Baldry I. K., et al. 2002, ApJ, 569, 582 
Becker R. H., White R. L., Helfand D. J., 1995, ApJ, 450, 

559 

Best P. N., Kauffmann G., Heckman T. M., Brinchmann 
J., Chariot S., Ivezic Z., White S. D. M., 2005, MNRAS, 
362, 25 

Best P. N., von der Linden A., Kauffmann G., Heckman 
T. M., Kaiser C. R., 2007, MNRAS, 379, 894 

Bruzual A. G., Chariot S., 1993, ApJ, 405, 538 

Cannon R., et al. 2006, MNRAS, 372, 425 

Condon J. J, Cotton W. D., Greisen E. W, Yin Q. F., 
Perley R. A., Taylor G. B., Broderick J. J., 1998, AJ, 115, 
1693 

Croom S. M., et al. 2002, MNRAS, 337, 275 
di Matteo P., Combes F., Melchior A.-L., Semelin B., 2007, 
A&A, 468, 61 

Dressier A., Oemler A., Poggianti B. M., Smail I., Trager 
S., Shectman S. A., Couch W., Ellis R., 2004, ApJ, 617, 
867 

Eisenstein D. J., et al. 2001, AJ, 122, 2267 
Eisenstein D. J., et al. 2003, ApJ, 585, 694 
Francis P. J., Hewett PC, Foltz C. B., Chaffee F. H., 

Weymann R. J., Morris S. L., 1991, ApJ, 373, 465 
Fukugita M., Ichikawa T., Gunn J. E., Doi M., Shimasaku 

K., Schneider D. P., 1996, AJ, 111, 1748 
Gebhardt K., et al. 2000, ApJ, 543, L5 p 
Gunn J. E., Stryker L. L., 1983, ApJS, 52, 121 
Hambly N. C, et al. 2001, MNRAS, 326, 1279 
Heckman T. M., Smith E. P., Baum S. A., van Breugel 



W. J. M., Miley G. K., Illingworth G. D., Bothun G. D., 

Balick B., 1986, ApJ, 311, 526 
Hine R. G., Longair M. S., 1979, MNRAS, 188, 111 
Johnston H. M., Hunstead R. W., Cotter G., Sadler E. M., 

2005, MNRAS, 356, 515 
Kriss G. A., 1994, in Crabtree D. R., Hanisch R. J., Barnes 

J., eds, Astronomical Data Analysis Software and Systems 

III Vol. 61 of ASP Conference Series, Fitting models to 

UV and optical spectra using SPECFIT in IRAF. ASP, 

San Francisco, pp 437-446 
McCarthy P. J., 1993, ARA&A, 31, 639 
Magorrian J., et al. 1998, AJ, 115, 2285 
Matthews T. A., Morgan W. W., Schmidt M., 1964, ApJ, 

140, 35 

Mauch T., Murphy T., Buttery H. J., Curran J., Hunstead 
R. W., Piestrzynski B., Robertson J. G., Sadler E. M., 
2003, MNRAS, 342, 1117 
Mauch T., Sadler E. M., 2007, MNRAS, 375, 931 
Owen F. N., Ledlow M. J., Keel W. C, 1995, AJ, 109, 14 
Rixon G. T., Wall J. V., Benn C. R., 1991, MNRAS, 251, 
243 

Roseboom I. O, et al. 2006, MNRAS, 373, 349 
Roseboom I. G. et al. 2007, The 2dF-SDSS LRG and QSO 

survey: the cosmic evolution of Luminous Red Galaxies, 

MNRAS, in prep. 
Sadler E. M., et al. 2007, MNRAS, 381, 211 
Sadler E. M., Jenkins C. R., Kotanyi C. G., 1989, MNRAS, 

240, 591 

Tadhunter C, Dickson R., Morganti R., Robinson T. O, 
Wills K., Villar-Martin M., Hughes M., 2002, MNRAS, 
330, 977 

Vanden Berk D. E., et al. 2001, ApJ, 122, 549 

Wills K. A., Tadhunter C. N., Robinson T. O, Morganti 

R., 2002, MNRAS, 333, 211 
York D. G., et al. 2000, AJ, 120, 1579 



© 0000 RAS, MNRAS 000, 000-000 



